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Abstract 

We consider proper colorings of planar graphs embedded in the 
annulus, such that vertices on one rim can take Q s colors, while all re- 



■ maining vertices can take Q colors. The corresponding chromatic poly- 

nomial is related to the partition function of a boundary loop model. 
Using results for the latter, the phase diagram of the coloring problem 
(with real Q and Q s ) is inferred, in the limits of two-dimensional or 
^ , quasi one-dimensional infinite graphs. We find in particular that the 

special role played by Beraha numbers Q = 4 cos 2 ^ for the usual chro- 
matic polynomial does not extend to the case Q ^ Q s . The agreement 
with (scarce) existing numerical results is perfect; further numerical 
checks are presented here. 



1 Introduction 

Let G = (V, E) be a planar graph embedded in the annulus. Let V B C V 
be the subset of vertices surrounding the face that contains the point at 
infinity. In other words, V s are the vertices on the outer rim of the annulus. 
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Place a spin variable o~i = 1,2, ... ,Q on each bulk vertex (i G V \V S ) and 
a boundary spin Uj = 1,2, . . . ,Q S on each boundary vertex (j G V s ). We 
suppose initially that Q s < Q, so that Q — Q s of the colors allowed for the 
bulk spins are forbidden for the boundary spins. 

The Potts model partition function Zg{Q, Q s ] v) — also known to graph 
theorists as the multivariate Tutte polynomial — can be defined through a 
slight generalization of the usual Fortuin-Kasteleyn expansion [T] 

(n \ ks(A) 

z G (Q,Q s ;v) = Q k{A) (77) Y[v e (1) 

ACE \ ) e£A 

where k(A) is the number of all connected components (clusters) in the graph 
induced by the edge subset A, and k s (A) is the number of connected compo- 
nents that contain at least one vertex from V s . In other words, Q s (resp. Q) 
is the weight of a cluster that contains at least one (resp. does not contain 
any) vertex in V s . The edge variables v = {v e } e( z E are related to the usual 
spin-spin couplings K e through the relation v e = exp(K e ) — 1. 

Note that in ([1]) there is no need for Q and Q s to be integers, nor do we 
have to impose the constraint Q s < Q. We shall henceforth promote ([!]) to 
the definition of the (boundary) Potts model [2]. 

In this paper we wish to study the problem of proper colorings of G, such 
that bulk vertices can have Q different colors, whereas boundary vertices can 
have only a subset of Q s colors. Adjacent vertices (of whatever type) are 
constrained to have different colors. The partition function Zg(Q,Q s ; —1), 
i.e. with all f e = —1, counting the number of such proper colorings is referred 
to as the boundary chromatic polynomial and denoted Pg{Q, Qs)- Note that 
Pg{Qi Q) is nothing else than the usual chromatic polynomial, which has 
been studied extensively in the literature [3]. 

We address in particular the issue of the phase diagram of Pg(Q, Qs) for 
"large graphs" — what is meant precisely by this will be discussed below. The 
main result is the location and nature of a series of phase transition (with 
corresponding behaviors of zeroes of Pq) occurring when one varies one or 
both of the parameters Q and Q s . We emphasize that most of our results are 
quite general, and do not depend on the detailed structure of the underlying 
graph G. 

In the usual case Q = Q s , one of the striking features of the chromatic 
polynomial is that for "large graphs" its real zeroes possess accumulation 
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Figure 1: For this graph (with the shaded vertices on the boundary) the 
chromatic polynomial continued from the region Q > Q s is Pq = {Q 2 — 3Q + 
3)Q s (Qs — 1)- It does not vanish for Q = nor Q = 1. Meanwhile the 
chromatic polynomial continued from the region Q < Q s is obtained — in this 
simple example — by exchanging Q and Q s in the above expression, and does 
vanish for Q = and Q — 1. 

points which belong to the magic set of Beraha numbers: 



Note that the first two such numbers are Q = and Q = 1, which are usually 
exact zeroes for finite graphs as well. One of the striking conclusions of our 
study is that the special role played by Beraha numbers is not very resistant 
to changing Q s . Depending on the problem one choses, there can indeed be 
accumulating zeroes at other special points of the real axis. 

It is important to realize that the definition of Pg{QiQs)i albeit very 
natural, can lead to counterintuitive features in particular when interpreted 
outside the initial domain of definition Q s < Q. For instance it turns out 
that for most graphs, Pg(Q, Qs) does not vanish when Q = or Q = 1, even 
though in that case there is no way — forgetting the boundary contribution — 
to color the bulk vertices with Q colors. The point is that in the original 
definition ([T|), spins belonging to clusters that touch the boundary receive 
a fugacity Q s , which initially is assumed smaller or equal to Q, but which, 
after continuation, can in fact be greater, hence "pumping" the number of 
colors in the bulk. Fig. [1] provides a simple example of this subtlety. 

One could define another chromatic polynomial starting from the situ- 
ation where Q < Q s . In terms of the subsequent cluster and loop model 
expansions, it would however be much less interesting. Indeed, in such a 
model, only clusters not containing any of the bulk spins would get the fu- 




for t = 2, 3, . . . . 



(2) 
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Figure 2: The transition system (shown as dashed lines) depends on whether 
a given edge e (show as a solid line) is present in [left panel] or absent from 
[right panel] the edge subset A C E in Q. 

gacity Q s , and thus only loops "glued to the boundary" would get a fugacity 
different from the bulk ones. This presumably would not affect the patterns 
of zeroes. 

The layout of the paper is as follows. In section [2] we relate the boundary 
chromatic polynomial to a loop model which was previously introduced in 
[2] and further studied in [1] . In section [3] the issue of the phase diagram is 
transposed into the setting of the Beraha-Kahane Weiss theorem [6] which 
we review. The necessary input for applying that theorem is supplied by 
an analytic continuation of the field theoretic results of [2], as explained in 
section HI Here we also arrive at the main results of the paper, which are 
the phase transition loci (I15p — (j!6p . All of this applies to the two-dimensional 
thermodynamic limit. However, the main results remain valid for quasi one- 
dimensional graphs, and we provide the necessary arguments in section [5j A 
few numerical validations of our results are given in section [6] after which we 
present our conclusions. 

2 Boundary loop model 

The cluster model ([1]) can obviously be defined for any graph G. However, 
when G is planar, the cluster model can be turned into a loop model on the 
medial graph M.. 

We recall that the medial (or surrounding) graph has a vertex standing 
on each edge e 6 E, and an edge between vertices standing on edges ei,e2, 
whenever ei, €2 are incident to a common vertex in V and surround a common 
face in G. 

A non-intersecting transition system on Ai is defined locally as in Fig. [2J 
Globally, this transition system is a set of loops — or cycles in the standard 
graph theoretical terminology — which separate clusters in G from their duals. 
By the existence of a point at infinity, the inside and outside of a loop are 




4 



well defined. A loop that contains at least one vertex of on its inside is 
called a boundary loop. A loop that is not a boundary loop is called a bulk 
loop. 

Let now £(A) be the total number of loops, and £ S {A) the number of 
boundary loops. By the Euler relation, one has k(A) = ^(£(A) + \V\ — \E\), 
so that 



Z G {Q,Q s -v) = Q\ y \l 2 Y.Q 

ACE 




where we have introduced x e = Q 1 ^ 2 v e . In other words, the weight of a bulk 
loop (resp. a boundary loop) is n (resp. n s ), subject to the relations 

Q = n 2 Q s = nn s (4) 

The boundary loop model ([3]) was introduced in [2J, and further studied 
in a more general setting in [I]. The emphasis in Refs. [2J H] was on the 
ferromagnetic case where all x e — 1. We shall see now that the generalization 
of these results to the antiferromagnetic region (with x e < 0) allows to infer 
the phase diagram of the boundary chromatic polynomial. 

3 Beraha-Kahane- Weiss theorem 

We wish to study the boundary chromatic polynomial in the thermodynamic 
limit where G becomes large (|V| — > oo). In general, one may take the limit 
\V\ — > oo through a recursive family of graphs Gn embedded in the annulus, 
of width W and circumference N, such that \V\ ~ NW and \V S \ ~ N. In 
particular one may think of strips of regular lattices (square, triangular,. . . ), 
but we emphasize that most of our results do not depend on the detailed 
structure of the graph, nor do they require that it be regular. 

In section 0] we take the width W oc N, so that the limiting graph G^ is 
two-dimensional, and the results [2] of conformal field theory (CFT) apply. In 
section [5], we consider instead W finite, so that Goo is quasi one-dimensional, 
and we shall see that the main results hold true in that case as well. 

In both cases one may think of the partition function Pg{Q, Q s ) as being 
built up by a transfer matrix, with time slices containing W spins. The 
structure of the transfer matrix has been discussed in details in |2j, and in 
particular it was shown that each of its eigenvalues Aj contributes to the 
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partition function with a non-trivial multiplicity Di that we shall refer to as 
an eigenvalue amplitude. Hence, 

P G (Q,Q S )=J2^) N ■ ( 5 ) 

i 

The fact that Di ^ 1 in general can be traced back to the non-local nature of 
the loops defining ([3]), and to the periodic boundary conditions in the time 
direction. 

We wish to study the phase diagram of the boundary chromatic polyno- 
mial by locating the accumulation points A of the partition function zeroes 
Pg{Qi Qs) — in the limit iV — > oo. Following Lee and Yang [3], this can be 
done by fixing one of the variables Q or Q s (or by imposing a fixed relation 
among Q and Q s ), and letting the remaining variable (henceforth denoted z) 
take complex values. 

Due to the form (jSJ) the Beraha-Kahane- Weiss (BKW) theorem [5] ap- 
plies. Let us call an eigenvalue Aj dominant at z if |Aj(z)| > |Afc(z)| for all k. 
The BKW theorem then states that (under very mild assumptions) 

• z G A is an isolated accumulation point iff there is a unique dominant 
eigenvalue Aj at z and the corresponding amplitude vanishes, a.i{z) = 0. 

• z G A forms part of a continuous curve of accumulation points iff there 
are at least two dominant eigenvalues at z. (In other words, z is the 
locus of a level crossing involving a dominant eigenvalue.) 

It is not in general clear to what extent CFT predictions apply to complex 
values of the parameters Q and Q s . But at least we can infer important 
information about the phase diagram by combining the BKW theorem [5] 
with the CFT results [2] for the special case of real parameter values. 



4 Phase diagram in the thermodynamic limit 

It is useful to parametrize the bulk and boundary loop weights as follows 

o ( \ sin((r + l)7re ) 

n = 2 cos(7re ) , n s = — (6) 

sm(r7reoJ 

defining the parameters eo and r G (0, — ). The continuum theory then has 
central charge 
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The range eo G [0, |) describes the usual ferromagnetic-paramagnetic transi- 
tion, corresponding to positive values of n and n s . 

We here need the analytic continuation into the range eo G (|, 1), where n 
and n s become negative. This range was referred to as the Berker-Kadanoff 
(BK) phase in [TJ. Inspecting Fig. [2] it is easy to see that (j3J) is invariant 
under a simultaneous sign change of n, n s , and x e . The BK phase therefore 
corresponds to negative values of x e , i.e., it describes a part of the antifer- 
romagnetic region of the Potts model. Its relevance to the chromatic line 
v e = —1 is due to the fact that the temperature variable v e is an irrelevant 
perturbation in the BK phase, in the sense of the renormalization group. 
The BK phase therefore controls, for any fixed Q G (0,4), a finite range of 
values v e . One may therefore expect that at least for Q < Q c , where Q c < 4 
is some lattice-dependent constant, the BK phase will control the chromatic 
line v e = —1. 

To give a little more substance to this general discussion, it is worthwhile 
to recall some exact information about the special cases of the square and 
triangular lattices. The standard Potts model (Q B = Q and v e = v) is then 
exactly solvable on the curves [8, 9J 

v 2 = Q (square lattice) 
v 3 + 3v 2 = Q (triangular lattice) (8) 

In view of the parametrization ([6]) it is more convenient to rewrite this as 

v = 2 cos(7re ) (square lattice) 

v = —1 + 2 cos ( (triangular lattice) (9) 

where e G [0, 1] for the square lattice and e G [0, |] for the triangular lattice. 

Both analytical and numerical studies of the Potts model with Q s = Q 
and either free or periodic transverse boundary conditions conclude that the 
critical exponents along the curves (ED for eo G [0, 1) are those predicted by 
the CFT. In particular, the central charge is ((71) as claimed. Moreover, the 
exponents for eo G (|, 1) are just the analytic continuations of those valid for 
the usual ferromagnetic regime eo G (0, |). This already strongly suggests 
that the critical properties for e G [0, 1) are lattice-independent (univer- 
sal) Q This conclusion is further corroborated by the so-called Coulomb gas 
approach |TD] to CFT. 

1 In the case of the triangular lattice, the range eo £ (1, §] describes a very different 
CFT 11 which we shall not need further in the present work. 
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Further studies have established that for each Q G (0, Q c ) the chromatic 
polynomial indeed renormalizes to the BK phase, with the following values 
of Q c for the square |T2j and triangular [13], [H] lattices 



Q c = 3 (square lattice) 

Q c = 3.8196717312- • ■ (triangular lattice) (10) 

We now return to the main objective of this section, which is to establish 
the critical behavior of the boundary chromatic polynomial. On general 
grounds, boundary conditions should not modify bulk RG flowsj^ Therefore, 
we expect the analytic continuation of the CFT results [2] to the range eo G 
(|, 1) to describe the critical behavior of the boundary chromatic polynomial 
for Q G (0,Q C ). 

It is convenient to set e = 1 — |, so that the BK phase corresponds to 
t > 2. The parameter r appearing in (jSj) is then constrained to r G (0, rzr)- 
We have 

n = —2 cos 



sin 



7T 
7 

(r(t-l)-l)ir 



™s = . V /.rJiwN 7 (11) 

sin 



^ r(t-l)^ 



In this parametrization, Q = n 2 [see Eq. (j3j)] is nothing else than the t'th 
Beraha number B t defined in (T5]). Real chromatic zeroes have long been 
known [3] to accumulate around B t for integer values of t > 2. One major 
motivation of this work is to show that the special role played by the Beraha 
numbers is destroyed by chosing Q s ^ Q. 

As explained in [2] the detailed transfer matrix structure implies that 
each eigenvalue appearing in (jSJ) is in fact an eigenvalue of a modified transfer 
matrix in which the number of loops winding around the periodic direction 
of the annulus (i.e., which are non-homotopic to a point) is fixed. Each 
eigenvalue can thus be labelled by the corresponding number of winding loops 



2 This is of course a subtle issue in cases such as this one, where the statistical models are 
not very physical. In fact, there are some boundary terms that can profoundly affect the 
behavior of flows in the Berker-Kadanoff phase — for instance those breaking the quantum 
group symmetry in the XXZ chain version of the models. For the boundary terms we 
are considering however — which can be described through the boundary Temperley Lieb 
algebra — no such "rogue" behavior seems to occur. 
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L — 0, 2, 4, . . ., as Aj . By the definition of the Potts model and the medial 
graph Ai, the corresponding number of winding clusters is L/2. In each 
sector with L > 0, the dominant eigenvalue corresponds to the outermost 
of the winding loops being constrained to be a boundary loop (i.e., we can 
restrict to what was called the "blobbed sector" in [21 H] ) . 

The existence of L winding loops corresponds in CFT to the insertion 
of a pair of so-called L-\eg operators Ol at the extremities of the strip; 
the extremities are subsequently glued together to form the annulus with 
periodic boundary conditions in the time direction. The asymptotic scaling 
for W ^> 1 of the dominant eigenvalue A^ in each sector L is then fixed by 
CFT as [15] 

A o L) ( * h A n ^ 

^ = + (12) 

where the dots on the right-hand side represent terms that decay faster than 
exp(— const /W). 

The constant Hl appearing in (fT2|) is the so-called conformal weight of the 
L-leg operator (in the "blobbed sector" ) whose value has been established in 
[21 [TH]. After the analytic continuation implied by the parametrization fllip . 
this reads 

h L = ± (L 2 - 2rL(t - 1) + (r 2 - l)(t - l) 2 ) (13) 
The corresponding eigenvalue amplitude has been derived rigorously in [2]: 

D - / 1 for L = , . 

L ~ \ n s U L ^(n/2) - U L _ 2 (n/2) for L > [ ' 

where Uj(x) is the j'th order Chebyshev polynomial of the second kind. 

Note that in the probabilistic regime (e G [0, |]) the continuum limit is 
dominated by L = 0. This is no longer true in the BK phase, where for any 
t there is at least one of the exponents hi taking negative values. The most 
negative exponent determines the most "probable" number of winding loops. 
This situation is clearly counter-intuitive from a probabilistic point of view, 
but it is made possible by the appearance of negative Boltzmann weights. 
Note also that the invariance of ([3]) under a simultaneous sign change of n, 
n s , and x e is not sufficient to make all weights positive. 

Using f[T2l , dominant level crossings of transfer matrix eigenvalues corre- 
spond asymptotically (for W ^> 1) to level crossings of the conformal weights 
hi. We can thus read directly from ( {TBI) the necessary and sufficient criterion 
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Figure 3: Conformal weights Hl as functions of the parameter r for the case 
t = 6. Dominant level crossings and vanishing dominant amplitudes are 
shown respectively as solid circles and crosses. 



for the second part of the BKW theorem. Indeed, level crossings involving 
the dominant L-leg sector occur when 

h L = h L+2 <£> r = — (15) 
with L < t — 1. In particular, Hl is the most negative exponent for r 6 

( L-l L+l \ 

\ t-i ' t-i r 

Similarly, the necessary and sufficient criterion for the first part of the 
BKW theorem is read off from ( fl4l) . Indeed, the amplitude of the dominant 
L-leg sector vanishes when 

D L = 0^r = j— j (16) 

with L = 2,4,6, .... 

These phenomena are illustrated in Fig. [3] for the case t = 6 (the Q = 3 
state Potts model). 

For any fixed n, phase transitions will therefore take place for r = s/(t — l) 
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and integer s £ (0, t\. The corresponding value of the boundary parameter is 



sin (f ) 



(17) 



For even s this corresponds to a vanishing dominant amplitude, and for odd 
s to a dominant level crossing. The corresponding value of the dominant 
exponent (|T3|) is h L = — for any even s, and /il = for any odd s. 

The iV — > oo limiting curve of accumulation points of partition function 
zeroes in the complex Q s plane (in the vicinity of the real Q s axis) can now 
be inferred from the BKW theorem: For even s one has an isolated real 
accumulation point, and for odd s a continuous curve of accumulation points 
intersects the real axis. 

In the example t — 6 of Fig. [31 the transitions at r = |,|,|,|,1,| 
correspond to the following numbers of boundary colors: Q s = 0, 1, |, 2, 3, oo. 

The discussion following ffTB]) has subsumed that we are interested in the 
phase diagram for fixed Q and varying Q s . But of course the criteria ([75]) - 
(|T6|) for phase transitions hold true for other situations as well. In particular, 
the following few useful cases correspond to simple relations between r and 
t: 

r = 1 



Q s 
Qs 
Qs 
Qs 
Q s 
Qs 
Qs 
Qs 



Q 
Q- 
Q- 

Q- 

o 

l 

2 

\Q 



r 
r 
r 
r 
r 
r 
r 



(t-2)/(t- 

(t-2)/(2f 

1/2 

l/(t-l) 
2/(t-l) 
(t + 2)/(2t 
t/(2t - 2) 



2) 



;is) 



For all of these, (Tl5!)-(fl6!) yield phase transitions located at integer values of 
t (i.e., at the Beraha numbers B t ), but this needs of course not be the case 
for more general choices of Q s . 



5 Quasi one-dimensional case 

We now turn to the quasi one- dimensional geometry where the circumference 
of the annulus iV — > oo, while is width W is kept fixed and finite. In that 
case, the possible number of winding loops is constrained by L < 2W. 
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Eq. 0141) for the eigenvalue amplitudes was in fact derived combinatorially 
for finite W, and so remains valid in this case. On the other hand, Eq. 0131) 
must be discarded, since its derivation supposed the validity of conformal 
field theory. However, the pleasant surprise is that even for finite W the 
dominant eigenvalues in the L and (L + 2) leg sectors cross exactly for the 
values of r and t given by OTol) . 

This coincidence follows from representation theory of the underlying 
boundary Temperley-Lieb algebra. While this algebra is semi-simple for 
generic values of the parameters, it admits families of degeneracy points 
where generically irreducible representations merge into larger indecompos- 
able representations. Results in [T7| guarantee that this occurs for finite 
values of W exactly at the same values that lead to the coincidences 0T5'|) of 
the conformal weights in the continuum limit. 

When r = 1 — in the original parametrization — this can be understood 
somewhat more easily by using quantum group representation theory [18J, 
as the generic U q (sl{2)) representations for sectors L and L + 2, of spin 
j = L/2 and j = L/2 + 1, merge into larger indecomposable representations. 
When r is integer larger than one, this can be explained similarly by the 
construction of section 5 in Ref . [2] . Indeed, there the effect of the boundary 
weight n s was obtained algebraically by adding r extra strands on the outside 
of the annulus, subject to the action of a certain symmetrizer. Thus, the 
boundary loop model ([3]) with r integer is a special case of the standard 
loop model in which only the weight n appears. The latter is known to have 
an U q (sl(2)) quantum group symmetry [18] . and this in fact implies that 
f fT5|) still holds true. The presence of exact coincidences at arbitrary r can 
maybe be interpreted in terms of some quantum group — the commutant of 
the boundary Temperley-Lieb algebra — but we will not discuss this here. 

The key results of section H] therefore remain valid, up to two subtle effects 
to be discussed below. 

To make this conclusion more accessible to readers unacquainted with 
quantum groups we turn to a numerical verification. Fig. H] shows the lead- 
ing free energies fa = — ^ log in the L-\eg sector, as functions of r in 
the parametrization QUI) , for four different values of W. The results were 
obtained for the square lattice in the diagonal geometry defined in [2] , along 
the curve Qj with e G (|, 1), i.e., within the BK phase. Results for other 
lattices would be similar, provided that one remains inside the domain of 
attraction of the BK phase. 

For each W, the dominant level crossings are seen to occur exactly as 
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Figure 4: Leading free energies fi in the sectors L = 0, 2, 4, 6 as functions of 
r G (0, |). The boundary loop model is here defined on the square lattice, 
along the BK critical curve, and the parameter t = 6. Four different system 
sizes {W = 8, 10, 12, 14) are shown, the largest size corresponding to the 
lowermost curves. The vertical lines are guides to the eye. 
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predicted by (IT51) . More generally, the r values singled out by (TI3I) (TTBT) 
are seen to be the loci of subdominant level crossings as well, as would be 
expected from an underlying quantum group symmetry. 

Fig. H] was made for the choice t — 6 (the Q = 3 state Potts model), so 
that it is the precise finite-size analogue of Fig. [3j Other, non-integer choices 
of t were found to lead to the same conclusions. 

We still need to discuss the two subtle effects referred to above. The 
first one is that if the annulus is too narrow (2W < \t\ ) to accomodate the 
number of legs required by dominant sector with the largest L predicted by 
(115]) . the corresponding level crossings will simply be absent, and the 2W-\eg 
sector will remain dominant for the corresponding values of the parameter r. 

The second effect is that Fig. H] gives clear evidence that when r becomes 
too large, there is an internal level crossing in each L-leg sector, visible as a 
cusp in the curves. To the right of these cusps the pattern of dominance may 
change. A detailed analysis of the loci of the cusps reveals that their position 
tends to r = 1 as W —>■ oo, independently of the value of L. Moreover, for 
r e (1, jzt) it is the L = sector that will be dominant for large enough W. 

6 Numerical verifications 

To conclude this paper, we wish to check that the predictions of sections HH5J 
agree with existing numerical results on the limiting curves A of chromatic 
zeroes. The goal of this comparison is furthermore to convince the reader 
that our results are: 

1. Lattice independent; 

2. Independent of v e , as long as we are in the BK phase; 

3. Correct for various choices of Q s . 

Fig. [5J shows the accumulation points A for the triangular-lattice chro- 
matic polynomial on an annulus of width W = 7. Transverse boundary 
conditions are free, so that Q s = Q. The agreement with the predictions 
( I15p - (|16p for the real accumulation points is perfect. There is one additional 
real accumulation point at Q C (W) = 3.4682618071 • ■ ■ which is a finite-size 
analogue of Q c discussed in section HI As W — > oo we expect Q C (W) — > Q c 
given by (ITU]) . 
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Zeros tri lattice L x = 7 R 
i 1 i 1 r 




i i i i i i i 

2 4 

Re(q) 

Figure 5: Zeroes in the complex Q plane of the triangular-lattice chromatic 
polynomial on an W x N annulus for W = 7 and N = 35, and their accu- 
mulation points as N — > oo. The boundary parameter Q s = Q. Taken from 
Figure 7 of Ref. ^9\. 
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Figure 6: Zeroes in the complex Q plane of a square- lattice Potts model 
along the curve (jSJ) on an W x iV annulus for W = 3 and N = 2Q, and 
their accumulation points as N — > oo. The boundary parameter Q s = Q — 1. 
Taken from Figure 5 of Ref . [20] . 
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Fig. E] shows the accumulation points A of partition function zeroes for 
a square-lattice Potts model along the curve (|SJ). The geometry is that of 
an annulus of width W = 3 with free transverse boundary conditions. How- 
ever, all vertices on the outer rim of the annulus are connected to an ex- 
tra exterior vertex. Therefore, the vertices on the outer rim (call them V s ) 
support spins which are effectively constrained to take only Q s = Q — 1 
different values (since they must be different from the value of the exterior 
spin). The partition function on the graph just described is therefore equal 
to QZg(Q, Q 8 — Q — 1; v e = ±\fQ) in our notation, where now G is just an 
ordinary annulus of width W, with no extra exterior vertex. 

Once again, the agreement with the predictions 0150 — (I16j) for the real 
accumulation points is perfect. In particular, it follows easily from the pre- 
dictions that the loci of isolated real accumulation points and curves of ac- 
cumulation points intersecting the real Q-axis are swapped between Figs. [5] 
and [6j Along the curve (jSJ) we would expect the BK phase to terminate only 
at Q c = 4. Thus, the phase transition corresponding to the largest possible 
L- sector becoming dominant is limited by the available width as L < 2W. 
This is again in perfect agreement with Fig. EJ Similar agreements are found 
with the numerical results for real accumulation points given in [20] in the 
case W = 4 (for which the complete limiting curve A was not computed). 

As a final check, we have computed the boundary chromatic polynomials 
with Q s = Q- 2onanH^x]V annulus for W = 2 and N = 100, for 
both the square and the triangular lattice. Their zeroes in the complex Q 
plane are shown in Fig. [3 The agreement between Eqs. f|T5|) - f|T6|) and the real 
accumulation points for the triangular lattice is striking. Notice in particular 
that we predict in general that only Beraha numbers of even order, viz. B t 
with t — 4, 6, 8, 10, . . ., can appear as accumulation points on the real Q axis. 
For the square lattice, the branch cutting the real axis at Q = 3 marks the 
termination of the BK phase, in agreement with (flQl) ; to the right of this 
branch one does not observe any further structure as expected. 

7 Conclusion 

To summarize, we have introduced a new graph coloring problem — the bound- 
ary chromatic polynomial — and identified the loci of phase transitions for real 
values of the parameters Q and Q s . Our results are lattice independent, and 
valid not only on the chromatic line but in the entire Berker-Kadanoff phase. 
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Figure 7: Zeroes in the complex Q plane of the Q s = Q — 2 boundary 
chromatic polynomials on an W x N annulus, with W = 2 and iV = 100, for 
both the square and the triangular lattice. The black vertical lines indicate 
the positions of the Beraha numbers (J2J). 
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While we have provided a number of striking numerical tests that validate 
our analytical predictions, we believe we have left ample space for further 
numerical investigations of the boundary chromatic zeroes for families of 
graphs embedded in the annulus. 

A straightforward extension of the work presented here would be to con- 
sider graphs on an annulus for which bulk spins can take values 1,2, ... ,Q, 
whereas spins on the outer (resp. inner) rim of the annulus are constrained 
to take values 1, 2, . . . , Q Q (resp. 1, 2, ... , Q\). Note that in the cluster ex- 
pansion analogous to ([I]), the number of spin values accessible to clusters 
touching both rims can be taken as a further independent variable Qb, not 
necessarily equal to min(Q ,Qi). 

Recent work on the corresponding two-boundary loop model furnishes 
the results for the eigenvalue amplitudes [I] and the critical exponents [21], 
analogous to ffl~3]) - ffT4"j) of this article. The phase diagram for real parameter 
values Q, Q Q , Q\, Qb can therefore be worked out along the lines presented 
here. 
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